Coulomb Gaps in One-Dimensional Spin-Polarized Electron Systems 
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We investigate the density of states (DOS) near the Fermi energy of one-dimensioiial spin- 
polarized electron systems in the quantum regime where the localization length is comparable to or 
larger than the inter-particle distance. The Wigner lattice gap of such a system, in the presence of 
weak disorder, can occur precisely at the Fermi energy, coinciding with the Coulomb gap in position. 
The interplay between the two is investigated by treating the long-range Coulomb interaction and 
the random disorder potential in a self-consistent Hartree-Fock approximation. The DOS near the 
Fermi energy is found to be well described by a power law whose exponent decreases with increasing 
disorder strength. 

PACS numbers: 73.20.Dx, 71.20.-b 
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The effects of electron interactions are particularly 
strong in one-dimensional (ID) systems, leading to Lut- 
tinger liquids or to various instabilities Random dis- 
order also has strong effects in ID systems: Regardless 
of disorder strength, all states are known to be local- 
ized with the localization length comparable to the mean 
free path The interplay of electron interactions and 
random disorder in low-dimensional systems is of great 
current interest H,^]. It has been shown recently that 
in the absence of disorder the long-range Coulomb inter- 
action between electrons reduces quantum fluctuations 
so that the ground state acquires quasi-long range or- 
der much close to a ID Wigner crystal 1^. This reduc- 
tion of quantum fluctuations is reflected in the density of 
states (DOS), which vanishes faster than any power as 
the Fermi energy /i is approached ||,^ : 

exp{-A[HE,/\E~fi\)]'/^}, 
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where A and Ec are appropriate constants. In the 
strongly localized regime, where the overlap of electron 
wave functions is negligible, electrons may be treated 
classically and the DOS exhibits a Coulomb gap of the 
form 0]: 



g{E) ^ [\n{Ej\E - iJi\)] 
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When the Thouless length is shorter than the localiza- 
tion length, the system is in the disordered Fermi liquid 
regime, where the overlap of electron wave functions can 
be significant. In this case the electron interactions can 
be included perturbatively within the weak localization 
theory js), and the flrst-order correction to the DOS near 
the Fermi level is given by H 



5g(E)^\E-^,\ 
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except in the limit E —i- fi, where the perturbation theory 
is expected to breakdown. 

In this work we employ a self-consistent Hartree-Fock 
(HF) method to investigate the DOS at zero tempera- 
ture in the quantum regime where the localization length 
is comparable to or larger than the inter-particle dis- 
tance . Both the long-range Coulomb interaction and 
disorder are expected to reduce quantum fluctuations: 
The Coulomb interaction pushes the system to the clas- 
sical limit, where quantum fluctuations can be neglected 
except for the redefinition of the strength of the impu- 
rity potential [ pdl , and disorder is expected to restore 
the Fermi liquid behavior jT^,^ . Thus a HF mean-field 
approximation may provide a reasonable description of 
the interplay between disorder and the Coulomb inter- 
action ||T^. This interplay is especially interesting with 
spin-polarized electrons since the Fermi wave vector is 
equal to 7r/a, where a is the period of the Wigner lat- 
tice. In the presence of weak disorder we expect that the 
processes leading to the Coulomb gap interact with the 
Bragg scattering leading to a gap in the DOS of a Wigner 
crystal. Such a spin-polarized system of electrons can be 
realized in organic chains or in quasi-one-dimensional 
quantum wires in a strong magnetic field ||l3|] . We flnd 
in the quantum regime that the DOS can be fitted well 
with a power-law, and that the exponent decreases as the 
strength of disorder is increased. 

A prototype model in which the interplay between dis- 
order and the Coulomb interaction can be investigated 
is a ID jellium model with N spin-polarized electrons, 
interacting with each other via the long-range Coulomb 
interaction in the presence of random impurities. We 
perform a finite-size calculation in a cell of length L with 



periodic boundary conditions. When the effects of the 
image charges are taken into account the Coulomb inter- 
action can be written as 



Vc{xi - X2) 
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where d is the cut-off length corresponding to the trans- 
verse dimension of the system [|| and e is the dielectric 
constant. For simplicity, we assume that the impurity at 
position Xi is characterized by the i5-function potential 
with random strength Wi. This gives the total impurity 
potential of Nj impurities in the form 
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where Wi and Xi are quenched random variables dis- 
tributed uniformly in the range [—Wmax/^, Wmax/'2] and 
in the interval [—L/2, L/2], respectively. The electron- 
electron interaction is treated self-consistently within a 
HF approximation in the momentum space. We expand 
the HF single-particle wavefunction for state a as follows: 



k 
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where the basis states are 



(t)k{x) = —j= exp{ikx) 

and the wavevector k takes integer multiples of 2n/L. 
The expansion coefficients cua satisfy the equations 



E 

fc' 



{k\HQ + Wi + U + VHF\k') Ck'a — faCfca, 
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where the matrix elements are evaluated by integrating 
the coordinate over the interval [—L/2, L/2]. Here Hq is 
the kinetic energy and the HF matrix elements are given 
by 



{k\VHF\k') ^Y.'^-^^^^°'\'^c\k' ,a) - {k,a\Vc\a,k') 
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where the first and second terms are the Hartree and ex- 
change potentials, respectively. The potential U is due to 
the uniform positive background charge and na is 1 /O for 
occupied/ unoccupied states a. The HF matrix elements 
in Eq. (||) depend on the expansion coefficients via 

{k,a\Vc\k',a) - 4^c.Ck,c.{k,k,\Vc\k\k2) (9) 
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{k,a\Vc\a,k') = cl^Ck,o.{k,k,\Vc\k2,k'), (10) 

ki ,k2 

where the Coulomb matrix elements are given by 

{kl,k2\Vc\k[k'2) = Sk,+k,M',+k',K{k, - k[) (11) 

with 
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In our numerical work, the electron number N is fixed 
to L/d, so that the inter-particle distance is equal to 
d. We measure the strength of the Coulomb interac- 
tion and impurity potential using dimensionless param- 
eters V = Vi/Ek and W = Vd/Ek, where Vi = e^/ed, 
Vd = Wmax/d^ and Ek = ?i^fco/2me with fco = 211 /d. 
The impurity number Ni is set equal to 5A^ while the di- 
mension of the Hamiltonian matrix is taken to be 2A^+1 
for W = 0.1, 3A^-|-1 for W = 0.4, and 57V-1-1 for W = 1.0. 
The convergence of the obtained results has been tested 
by increasing the dimension of the matrix. Figure 1 dis- 
plays the probability distribution function of an electron 
at the Fermi energy in strongly localized [W = 1.0), 
intermediate {W = 0.1), and disorder free {W = 0) 
regimes. We see that the localization length decreases 
with the disorder strength, and particularly, in the ab- 
sence of impurities, observe ID crystalline order rather 
than an electron liquid. This mean- field result is in quali- 
tative agreement with the result obtained through the use 
of the bosonization technique, which has predicted quasi- 
long-range-order with the correlation function decaying 
slower than any power, regardless of the strength of the 
long-range Coulomb interaction [||. At W — 0.1 the 
electrons are in the quantum regime since the localiza- 
tion length extends over several times the inter-particle 
distance. We note that even at W = 1.0 the electrons 
the localization length is comparable to the inter-particle 
distance. 

The average DOS integrated over the system size is 
given by 
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and 



where AD{E) represents the number of eigenenergies in 
the energy range of width AE around energy E and 
the double angular brackets denote the disorder aver- 
age. Here we set the value of AE to be 0.05 in units 
of Ek, and compute the DOS by means of Eq. (p^. Fig- 
ure 2 displays g{E) at = 0.1 for five different system 
sizes L/d = 30, 40, 50, 60, and 80. We observe that a 
large reduction of the DOS occurs near the Fermi energy. 
Since there is an energy range where no significant DOS 
is present, we need to average over many disorder realiza- 
tions to get accurate values of the DOS. We find that for 
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W = 0.1 the number Nu of disorder realizations between 
1001 and 3844 is quite sufficient. To analyze the shape 
of the gap quantitatively, we fit the DOS near the Fermi 
energy to a power law: 



giE)^\E-^l\■ 
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where the DOS g{E) is related to the obtained (average) 
DOS g{E) via 



m = 



1 

Xe 
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This fit of the DOS for W = 0.1 is shown in Fig. 3. 
The power law fits our numerical results quite well, with 
/i = -0.29 ± 0.02 and a = 5.80 ± 0.16. For compari- 
son, we have also fitted g{E) to a logarithm, given by 
Eq. (||), only to find substantial disagreement: The func- 
tion, Eq. (|^), is a sublincar function of E while the nu- 
merical data form a superlincar function; see Fig. 2. 

Figure 4 shows the system size dependence of the ex- 
ponent a for two different values of disorder strength, 
W = 0.1 and 0.4. The values of {L/d, Nd) for W = 0.1 
are as in Fig. 2, while for W = 0.4 four new values 
are chosen. The numbers of disorder realizations for 
W = 0.4, which have been chosen between 199 and 
477, can be much smaller than those for T4^ = 0.1 since 
there are many more states near the Fermi energy due 
to stronger disorder. The obtained results of the size de- 
pendence tend to bend up in the large system, and yield 
finite extrapolated values of a in the limit L/d — s- oo, 
which provides an evidence for the algebraic gap forma- 
tion in the thermodynamic limit. As expected, the expo- 
nent a decreases with W, implying that the gap softens 
as the strength of disorder is increased. On the other 
hand, we expect that in the very clean limit {W <C 0.1) 
the power law will fail, leading to the gap of the form in 
Eq. (I). 

In summary, we have investigated the density of states 
of interacting spinless electrons in the presence of random 
disorder in one dimension. We find numerical evidence 
for the formation of a gap near the Fermi level: The 
density of states can be fitted with a power law in the 
quantum regime, where the localization length is compa- 
rable to or larger than the inter-particle distance. The 
magnitude of the exponent describing the algebraic gap 
decreases with increasing disorder. 
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FIG. 1. Probability distribution function, in units of d ^, of an electron at the Fermi energy for three different values of 
disorder strength. 
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FIG. 2. Density of states as a function of energy for various system sizes. The data marked by diamonds(O), plus signs(+), 
squares(n), crosses(x), and triangles(A) correspond to (L/d, Nd)={30, 1924), (40, 3740), (50, 3108), (60, 3844) and (80, 1001), 
respectively. The disorder strength W and the electron density n are fixed to 0.1 and l/d, respectively. 



FIG. 3. Power-law fit of g{E) in a system of size L/d = 60 and disorder strength W = 0.1. 



FIG. 4. Exponent a as a function of the inverse of the system size for two values of disorder strength, W = 0.1 and 0.4. 
While the parameters {L/d, No) for W = 0.1 are the same as those in Fig. 2, for W = 0.4 we have used {L/d, Nd)={30, 477), 
(40, 296), (50, 252), and (60, 199). The electron density is fixed to l/d. The lines are guides to the eye. 
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